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Abstract. Galactic cosmic rays are believed to be generated 
by diffusive shock acceleration processes in Supernova Rem- 
nants, and the arrival direction is likely determined by the 
distribution of their sources throughout the Galaxy, in par- 
ticular by the nearest and youngest ones. Transport to Earth 
through the interstellar medium is expected to affect the cos- 
mic ray properties as well. However, the observed anisotropy 
of TeV cosmic rays and its energy dependence cannot be ex- 
plained with diffusion models of particle propagation in the 
Galaxy. Within a distance of a few parsec, diffusion regime 
is not valid and particles with energy below about 100 TeV 
must be influenced by the heliosphere and its elongated tail. 
The observation of a highly significant localized excess re- 
gion of cosmic rays from the apparent direction of the down- 
stream interstellar flow at 1-10 TeV energies might provide 
the first experimental evidence that the heliotail can affect 
the transport of energetic particles. In particular, TeV cosmic 
rays propagating through the heliotail interact with the 100- 
300 AU wide magnetic field polarity domains generated by 
the 1 1 year cycles. Since the strength of non-linear convec- 
tive processes is expected to be larger than viscous damping, 
the plasma in the heliotail is turbulent. Where magnetic field 
domains converge on each other due to solar wind gradient, 
stochastic magnetic reconnection likely occurs. Such pro- 
cesses may be efficient enough to re-accelerate a fraction of 
TeV particles as long as scattering processes are not strong . 
Therefore the fractional excess of TeV cosmic rays from the 
narrow region toward the heliotail direction traces sightlines 
with the lowest smearing scattering effects, that can also ex- 
plain the observation of a harder than average energy spec- 
trum. 



1 Introduction 

During the last decades, galactic cosmic rays have been 
found to have a small but measurable energy dependent 
uneven arrival direction distribution, with a relative ampli- 
tude of order 10~ 4 — 10~ 3 . This anisotropy was observed 
in the northern hemisphere from energie s of tens to sev- 
eral hundreds GeV wi t h muo n detectors dNagashima et all 
" 20101) . and in the multi-TeV en- 



1998; Munakataetal 



ergy ra nge with Tibet AS7 a rray ( Amenomori et al.l 2006 



201 lab. Super-Kamioka nde dGuillian et all 120071) Mila 



gro dAbdoetajl l2009h and ARGO-YBJ dZhangl 120091: 
IShuwangu201 ll) . An anisotropy was also observed at an en- 
ergy i n excess of about 10 TeV with the E AS -TOP shower 
array (Ag lietta et al.l 120091) . Recently similar observations 
were reported in the southern hemisphere at en ergies of IP's 
to 100's TeV with the IceCube Observatory (Abbas i et al.l 
20101 201 lab . While at ~ 10 TeV the anisotropy appears to 
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be topologically connected to the GeV- TeV observations in 
the north, above about 100 TeV the global anisotropy per- 
sists b ut with a different p hase, consistently with the results 
from lAglietta et al.l (12009). The top panel of Fig.[T]shows the 
combined map in equatorial coordinates of relative intensity 
of cosmic ray arrival direction distribution observed by Ti- 
bet AS7 at about 5 TeV in the northern hemisphere, and by 
IceCube at about 20 TeV in the southern hemisphere. 

The arrival distribution of sub-TeV cosmic rays re- 
vealed th e existence of tw o kinds of anisotropies, as dis- 
cussed in iNagashima et al. dl998l) . One is a global non- 
dipolar anisotropy with relative excess approximately cen- 
tered around equatorial right ascension of hr, and an in- 
creasing amplitude up to 1-10 TeV energies. The other is 
a directional excess region confined in a cone of half open- 
ing angle of 68° from right ascension of about 6 hr, and ob- 
served for energies below TeV. This region covers a portion 
of the sky that includes the direction of the heliospheric tail 
(or heliotail), which is the region of the heliosphere down- 
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was r eported by IceCube, in the southern hemisphere, as 



well dAbbasi et all 1201 lbb . The bottom panel of Fig. [T] 



Milagro + IceCube TeV Cosmic Ray Data (10" Smoothing) 




significance 



Fig. 1. Top panel: Map in equatorial coordinates of the 
relative intensity of the cosmic ray arrival distribution as ob- 
served by th e Tibet AS7 at about 5 TeV in the northern hemi- 
sphere (from lAmenomori et al. I fcOllab ) and by the IceCube Ob- 
servatory at a median energy of 20 TeV in the southern hemi- 
sphere (from IXbbasi et all d2010l) ). Bottom panel: Map in equa- 
torial coordinates of the statistical significance of the cosmic ray 
arrival direction distribution as obs erved by Milagro a t about 1 TeV 
in the northern hemisphere (from lAbdo et ail d2008l) ) and by the 
IceCube Observat ory at a median energ y of 20 TeV in the southern 
hemisphere (from lAbbasi et all d201 lbh ). In this map features with 
angular extension larger than 30° — 60° are filtered out. 



stream the interstellar wind delimited within the heliopause, 
i.e. the bo undary that separates the solar win d and interstellar 
plasmas dlzmodenov and Kallenbachlkood) . Its origin was 
therefore attributed to some unidentified anisotropic process 
occurring in the heliotail, and thus it was called "tail-in" ex- 
cess. 

At the higher TeV energies, while the tail-in broad ex- 
cess becomes sub-dominant, the global anisotropy shows 
evidence of statistically significant small angular structures 
from the same direction in the sky. In particular, using ex- 
perimental techniques in the attempt to isolate relatively lo- 
calized excess or deficit regions of events that overlap over 
the smooth global anisotropy, angular features of order 20°- 
30° were discovered. Two separate highly significant lo- 
calized fractional excess regions of cosmic ra ys were re- 
porte d in the northern hemisp here by Milagro (|Abdo et all 
20081) . and al so by Tibet AS - y dAmenomori et all|2007l) and 
ARGO-YBJ dVernetto et all l2009t lluppal l201ll) . The ob- 
servation of a small scale anisotropy at multi-TeV energies 



shows the combined map in equatorial coordinates of statis- 
tical significance of the cosmic ray arrival direction distribu- 
tion where only features with angular extension smaller than 
about 60° are visible. Such small scale features lay in the 
same portion of the sky where the tail-in excess was domi- 
nant at lower energy, especially the one toward the heliotail 
direction with equatorial coordinates {a, 5) « (5 hr, +17°). 

At an energy in excess of about 100 TeV, where the 
anisotropy has a different topology than at lower energy, cos- 
mic ray particles are hardly influenced by the heliosphere and 
its elongated tail, and their arrival direction might hold in- 
formation on the Local Interstellar Magnetic Field (LIMF) 
on a larger scale. If the extended heliotail induces a signif- 
icant perturbation in the local interstellar medium, that can 
affect the arrival direction of multi-TeV cosmic ray parti- 
cles, then the anisotropy can be considered as an indirect 
probe of how t he LIMF influences the he liospheric bound- 
ary itself (see |De^mtiandLazarianj (|201l|)). Moreover, cos- 
mic rays below about 10 TeV are expected to be influenced 
by magnetic fields inside the heliotail as well. The con- 
current effects of magnetic reconnection and scattering pro- 
cesses might be able to explain some observations, although 
more experimental results and further developments in he- 
liospheric Magneto Hydro-Dynamic (MHD) simulations are 
needed for better constrain models. 

The origin of the cosmic ray anisotropy, its persistence 
in a wide energy range and its angular structure, is cur- 
rently subject of debate. In this paper we will briefly report 
the interpretations provided by various authors (in Sec. |2), 
with an emphasis on a possible phenomenological connec- 
tion between the broad tail-in excess of sub-TeV cosmic 
rays and the localized fractional excess of multi-TeV cos- 
mic rays from the direction of the heliotail. We'll then de- 
scribe the magnetic field structure in the heliotail as shaped 
by solar cycles and rotation in Sec. [3] An overview on mag- 
netic reconnection processes is given in Sec. |4] with an em- 
phasis on stochastic reconnection, assumed to contribute to 
the origin of the anomalies observed toward the heliotail. 
Sec. [5] addresses the proposed mechanism of cosmic ray re- 
acceleration in the heliotail and its effective influence in re- 
lation to scattering processes. 



2 Cosmic Rays Anisotropy 

The origin of cosmic ray anisotropy is still unknown. The rel- 
ative motion of the solar system with respect to the cosmic 
ray plasma rest frame (for instance due to galactic rotation) 
would p roduce a dipolar aniso t ropy in the direction of the 
motio n dCompton and Gettinpl 1 1935c iGleeson and Axford , 
1968). Such Compton-Getting effect was not singled out 
from observations yet, inducing to a possible conclusion 
that the bulk of galactic cosmic rays co-rotates with the so- 
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lar system (lAmenomori et all 120061) . Moreover, in the sce- 
nario where galactic cosmic rays are accelerated in super- 
nova remnants, their arrival direction should have a rela- 
tive excess toward the galactic center, i.e. the line of sight 
with the larger expected number of sources. On the other 
hand the nearest galactic sources would dominate the ob- 
served arrival distribution, and changes in anisotropy ampli- 
tude and phase with cosmic ray energy can arise as a nat- 
ural consequence of the stoc hastic nature of their sources in 
the local interstel l ar me dium (Erlv krnand Wolfendalel 2006 : 
Blasi and Amatol l2012h . 

Propagation properties of cosmic rays in the interstellar 
medium are likely to have an importan t role in shaping the 
anisotropy as well dBattaner et all 12009). For instance, a sce- 
nario where the large scale anisotropy is linked to diffusion 
of cosmic rays through the LIMF connecting the solar sys- 
tem to the interstellar medium outside the local interstellar 
cloud (wh ere the solar system currently r esides), was pro- 
posed by lAmenomori et all (120071 1201 lbb . This model ac- 
counts for the apparent quadrupol ar contributio n observed 
with the large scale anisotropy. In iFrischl (l201lh it is noted 
that the tail-in excess region, besides including the heliotail 
direction, is centered around the direction of the LIMF, there- 
fore linking its origin to their propagation deep inside the tail 
or to streaming along the LIMF or the S 1 sub-shell of Loop 
I superbubble. 

Within a distance of a few times the mean free path, dif- 
fusion regime breaks down and propagation of cosmic rays 
depends on their interaction with the turbulence ripples of 
the LIMF. Even though observations suggest that the LIMF 
is coherent over scales of about 100 pc, they also imply vari- 
ations in field directions of less than 3 0° - 40°, that can be 
attributed to turbulence (Frisch, 201 1 ). Scattering of TeV- 
PeV cosmic ray particles with the turbulent interstellar mag- 
netic field within the mean free path (i.e. a few 10's pc) can 
generate intermediate and small sc ale perturbations over a n 
underlying large scale anisotropy dGiacinti and Sigll l201lh . 
The observed anisotropy structure, therefore, could be used 
to infer turbulence properties of the LIMF. At energies be- 
low about 100 TeV the proton gyro-radius is a few thou- 
sands AU, thought to be comparable to t he length of the he- 
liotail (llzmodenov and Alexashovl 120031) . At these energies 
cosmic ray anisotropy is likely influenced by the extended 
and turbulent heliospheric magnetic field, and localized fea- 
tures in arrival direction can arise from the scattering of en- 
ergetic cosmic ray particles with the h eliospheric magnetic 



field ordered by the LIMF direction (Desi ati and Lazarian 
1201 lh . 

Another model aimed to explain the origin of the TeV 
small scale anisotropic features, appeals to the observation 
that the two localized excess regions in the northern hemi- 
sphere are seemingly close to the so-called Hydrogen De- 
flection Plane (HDP), which is the plane containing the di- 
rections of the interstellar flow and of the magnetic field 
upstream the heliospheric nose (lAmenomori et all 1201 lbb . 



According to this model cosmic rays propagating along the 
heliotail within the HDP are bent by the heliospheric mag- 
netic field so that two localized excess regions are formed 
symmetrically separated with respect to the direction of the 
heliotail on the HDP. This implies that the heliospheric mag- 
netic field between about 70 AU and 340 AU along the he- 
liotail is responsible for the two localized regions observed 
in the northern hemisphere in the energy range between 4 
and 30 TeV. The heliospheric magnetic field has a com- 
plex structure determined by the combined effects of the 26 
day rotation period of the S un and of the 1 1 year solar cy- 
cle dPogorelov et al. . 2009a ). This complex time-dependent 
magnetic field structure should produce an observable time 
variability in the relative intensity and position of the local- 
ized fractional regions over an 1 1 year period. 

Some other models rely on an astrophysi c al ori - 
gin 



of the observation. In Salvi 



Drurv and A haroniar] d2008h : ISalvatil d2010l) it is noted that 



tiandSacco (2008); 



the two observed localized excess regions in the northern 
hemisphere, surround the present day apparent location of 
Geminga pulsar. The supernova that gave birth to the pulsar 
exploded about 340,000 years ago, and the accelerated cos- 
mic rays might have propagated along interstellar magnetic 
fields connecting the region of Geminga to Earth. Since noth- 
ing or very little is known of the local interstellar medium 
properties, cosmic ray diffusion is not sufficiently constraint 
to provide a coherent scenario that can explain the observa- 
tions without considerable fine tuning. 

Due to the coincidence of the most significant localized 
excess observed by Milagro with the heliotail, it is possi- 
ble that we are seeing the effects of neutron production in 
the gravitational ly focussed tail of the inter stellar material, 
as suggested by iDrury and Aharonian d2008l) . Cosmic rays 
propagating through the direction of the tail interact with 
matter and magnetic fields to produce neutrons and hence 
a localized excess of cosmic ray in that direction. But while 
the target size has about the right size compared to the de- 
cay length of multi-TeV neutrons (~ 0.1 pc), the increase of 
the gravitating matter density is too low to account for the 
obse rved excess. 

In Malkov et al. (2010) it is proposed that cosmic rays 
emitted by a source (like a supernova remnant for instance) 
within a few 100 pc, are scattered by a strongly anisotropic 
Alfven wave spectrum formed by the turbulent cascade 
across the local field direction. Cosmic rays with small pitch 
angle with respect to the local interstellar magnetic field un- 
dergo the highest scattering, thus producing a faint localized 
excess region. An outer scale of the interstellar medium tur- 
bulence of about 1 pc would explain the observations. 

The fractional excess relative to the cosmic ray back- 
ground observed by Milagro in the direction of the heliotail 
is ~ 6 x 10~ 4 , i.e. about 1/10 the amplitude of the 
global anisotropy at TeV energy. This is comparable to the 
amplitude that the broad tail-in excess woul d have if extrap- 
olated from 100's GeV to TeV energies dNagashima et al 
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Fig. 2. A schematic meridional view of the sectored heliospheric 
magnetic field arising from the 11 -year solar cycle polarity rever- 
sals. The subsonic solar wind pushes the sectors across the termi- 
nation shock along the heliotail, compressing them to each other 
and toward the heliopause. The turbulence that is expected to per- 
turb the heliospheric magnetized plasma, determines the thickness 
of the current sheet outflow regions. While their length depends on 
the average magnetic field geome t ry and by the turbulence l evel. 
Adapted from lNernev et al.ldl99^ : lLazarian and Opherl(l2009|) . 



1998). Such an excess was found to be consistent with 
hadronic cosmic rays with an energy spectrum of the form 
N(E) ~ E~~<e~ E/E <= with spectral index 7 <2.7 (i.e. flatter 
than the average cosmic ray spectrum) at 4.6 a le v el, an d 
a cut-off energy E c = 3 - 25 TeV dAbdo et all |2008|) . 
A simil ar spectral hardening was observed by A RGO- 
YBJ dDi Sciascio and for the ARGO-YBJ Collaboration , 



20121) . 

In this paper we discuss the scenario where the excess re- 
gion of cosmic rays from the direction of the heliotail ob- 
served from a few tens GeV to about 10 TeV, is generated by 
re-acceleration processes of a fraction of energetic particles 
propagating through magnetic reconnection regions along 
the heliotail. A concurrent contribution from scattering with 
the turbulence ripples of the heliospheric magnetic field can- 
not be excluded, especially in relation to the other observed 
localized fractional excess re gions, although this p ossibility 



is the topic of another paper (IDesiati and Lazariani 1201 II) 



3 Magnetic field structure at the heliotail 

The motion of solar system through the local partially ion- 
ized medium^ produces a comet-like interface due to the 
solar wind plasma advected downstream by the interstellar 
flow, called the heliosphere. A termination shock, where the 
solar wind pressure equals that from the interstellar flow, is 
formed at approximately 100 AU from the Sun. The interface 
separating interplanetary and interstellar magnetic fields, 
called heliopause, is at a distance of approximately 200 AU 



the solar system is located at the edge of the so-called local 
interstellar cloud, which is part of a complex cloudlet structure ex- 
pandin g from the Scorpion-Centaurus Association (see lFrisch et al.l 



in the upstream dire ction, and it may extend downstrea m sev- 
eral thousands AU (llzmodenov and Kallenbach, 2006) where 
it could be about 600 AU wide dPogorelov et al.[|2009al) . The 
LIMF drapes around the heliosphere impri nting a deforma 



tion th at affects its internal structure as well (Pog orelov et al 
2009b). The heliospheric magnetic field has been studied 
with detailed MHD simulations, where the effects from the 
26 day solar rotation and the 1 1 year solar cycle were con 



sidered dPogorelov et all l2009al) (see also IScherer and Fahr 
(2003)). Over solar cycles the magnetic field polarity is re- 



versed every 1 1 years, generating unipolar regions dragged 
along the heliotail by the ~ 100 km/s solar wind jParker, 
In particular these magnetic regions grow to their 



197 



maximum latitudinal extent during solar minimum (about 
200-300 AU in size) and reduce to zero at solar maximum, 
when the heliospheric plasma is dominated by the strongly 
mixed polarity domains (about .1-1 AU in size) from so- 



lar rotation dNernev et all Il995l) . Due to the tilt of the 
solar magnetic axis with respect to its rotation axis, the 
unipolar regions are thinner at lower latitudes (as shown in 
Fig. [2). Therefore the tailward line of view is dominated 
by the finely alternating magnetic field, while along sigh- 
lines away from it the magnetic domains have larger size. 
MHD numerical simulations show that the sectored unipo- 
lar magnetic field regions can propagate for several solar cy- 
cles before they dissipate into the local interstellar medium. 
The corresponding periodic variations on the heliospheric 
plasma induce changes in the magnitude of the Alfven ve- 
locity by about 20%, and o f the magnetic field by about 
25% dPogorelov et alll2009ab . 

There is observational evidence that the plasma 
in t he heliosheath has Reyno lds number Re ~ 10 14 
(see lLazarian and Opherl (120091) and references therein), 
meaning that the strength of non-linear convective processes 
at the largest scale is more important than the damping 
viscous processes in the dynamics of the flow. We expect a 
similarly high Reynolds number in the inner heliotail as well. 
In such conditions it is very unlikely that plasma flow stays 
laminar, and the downstream motion in the heliotail is likely 
turbulent. In addition, the presence of neutral atoms in the 
partially ionized local cloud medium (where the solar system 
is moving) is essential for the dynamics of the heliosphere 
and LIMF interaction. Charge-exchange processes between 
the interstellar inflowing neutral atoms and the outflowing 
solar wind protons can produce Rayleigh-Taylor type insta- 
bilities on the heliopause with amplitude of a few tens AU 
and o ver a time scale of a few hundreds years dLiewer et al. , 
Il996h . Also in a model of plasma-neutral fluid coupled via 
collision and charge-exchange processes, it is found that 
such non-linear coupling leads to alternate growing and 
damping of Alfvenic, fast and slow turbulence modes, at 
L ^100's AU scale and with evolution tim e longer than 
inertial time L/V A dShaikh and Zankl l20ld) . with V A the 
Alfven velocity. Such modulations can propagate on the 
heliopause, producing ripples along the heliotail that can 
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penetrate deep inside the heliosheath and propagate outward 
into the local interstellar medium. Therefore, although 
more investigations are needed in order to understand the 
detailed plasma properties in the heliotail and its outer 
boundary, it is reasonable to assume here that magnetic 
fields in the heliotail are weakly stochastic, and likely 
reconnecting as the gradient in solar wind advects magnetic 
field lines closer to each other. The Alfven velocity of the 
turbulence in the heliotail is expected to be approximately 
40-70 km/s, with the actual value dependin g on the locatio n 




Sweet-Parker model 



within the sectored magnetic field (Pogorelo v et al. , 2009a). 



This is smaller than the solar wind speed downstream the 
termination shock, therefore magnetic reconnection in the 
heliotail is not expected to change the overall magnetic field 
structure. Nevertheless, the effects of turbulence are very 
important from the point of view of magnetic reconnection 
and the particle acceleration that it entails. 

Simulations of the magnetic fields in the heliotail are ex- 
tremely challenging due to its extension and to the com- 
plex interaction with the interstellar wind and between he- 
liospheric magnetic field and the LIMF, but mainly because 
there is currently no direct data collection from this remote 
portion of the heliosphere. Future refinements of MHD sim- 
ulations will provide higher resolution mapping of the helio- 
tail and of the plasma properties, that will help improving our 
knowledge of its effects on TeV cosmic ray propagation. 



4 Stochastic magnetic reconnection 

Astroph ysical p l asmas are often highly ionized and mag- 
netized ( Parkerl 1970l) . and they undergo dissipative pro- 
cesses, which annihilate the magnetic fields and convert elec- 
tromagnetic energy into plasma energy. Due to these pro- 
cesses, plasma from regions of a given polarity becomes 
magnetically connected to that of opposite polarity: this is 
when magnetic reconnection occurs. However, reconnection 
speed, and therefore the rate at which magnetic energy is 
converted into plasma energy, is too small to be important 
for acceleration of energetic particles, unless the effects of 

plasma resistivity are negligible. 

In the Sweet-Parker model of reconnection ( Sweeti 1958 



IParkeil 1 1957b the outflow is limited within the transition 
zone A, which is determined by Ohmic diffusivity (see top 
of Fig. [3j. In this model reconnection speed is smaller 
than the Alfven velocity of the plasma by a factor equal to 
S^ 1 / 2 = (LVa/i])^ 1 ^ 2 , where S is the Lundquist number, L 
the length of the current sheet, Va the Alfven speed and r\ is 
the Ohmic resistivity of the plasma. The length of the current 
sheet is determined by the extent of magnetic flux tubes that 
get in contact. Although the properties and dimensions of the 
heliotail are not well constraint yet, it is possible to state that 
the extension of current sheets between sectored heliospheric 
magnetic field in t he heliotail coul d lay be tween about 100 
AU and 300 AU ( Pogorelov et all l2009al) . Assuming the 




Turbulent model 



blow up 



Fig. 3. Upper plot: Sweet-Parker model of reconnection JSweetl 
ll958l : |Parkerlll957l) . The outflow region is limited within a thin 
transition zone A between the reversed magnetic field lines, which 
depends on plasma resistivity. The other scale is an astrophysi- 
cal scale L> A. Middle plot: Reconnection of weakly stochastic 
magnetic field according to lLazarian an d Vishniacj| |l999l) . The out- 
flow region is determined by the diffusion of magnetic field lines, 
which depends on the field stochasticity. Lower plot: An individ- 
ual small scale reconnection region. The reconnection over small 
patches of magnetic field determines the local reconnection speed. 
The global reconnection speed is substantially larger as many in- 
dependent patches come together. The bottleneck for the process 
is given by magnetic field wandering and it gets comparable to L 
as th e turbulence injection velocity approaches the Alfvenic one. 
From lLazarian et al.ld2004l) . 



same plasma properties as in the helios heath closer to the he- 
liosph eric nose, S is about 10 12 -10 13 dLazarian and Ophet , 
2009). Therefore, the corresponding reconnection speed 
for the Sweet-Parker model is several orders of magnitude 
smaller than the Alfven velocity. In fact in this case plasma 
collected over the size L should be ejected with speed ~ Va 



through the outflow region of thickness A = LS^ 1 ^ 2 , i.e. 
much smaller than the length of the current sheet. It is 
the large difference between L and A that makes reconnec- 
tion slow and unlikely to produce any effect on the plasma. 
The major consequence of such a model is that reconnection 
speed is limited by Ohmic resistivity of the plasma (see right 
panel of Fig . |4j» . Since most astrophysical plasmas have very 
low resistivity, a magnetic reconnection mechanism such this 
would not produce observable effects. 

On the other hand, various observations suggest that re- 
connection, when it occurs, can be fast in some circum- 
stances. For instance the development of solar flares sug- 
gests that magnetic reconnection should be initially slow in 
order to ensure the accumulation of magnetic flux, and then 
suddenly becomes fast in order to explain the observed fast 
release of energy. Fast reconnection would require L ~ A, 
meaning that the region over which magnetic flux tubes in- 
tersect is comparable to the size of the outflow region. This 
can be achieved by increasing the outflow region beyond 
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Fig. 4. Left panel: dependence of the reconnection speed Vr on injection power Pi n j. Rig ht panel; dependence of the reconnection speed 
Vr on the uniform resistivity rj lt . Open symbols are for Sweet-Parker reconne ction scenario (IS weetll 1 958l; IParkeJ, 1 1 957t) . and filled symbols 
are for weakly stochastic reconnection scenario (ILazarian and VishniacL[T999h . See lKowal et alj l2009h :l Lazarian et alj d201 ll) . 



the prediction of Sweet-Parker model, or by making L as 
small as the Ohmic diffusion region, so that magnetic field 
lines re connect in an "X -point". In this X-point collisionless 
model (Petschek, 1964) . reconnection speed does not depend 
on the resistivity and it is of the order of Alfven velocity of 
the plasma. On the other hand X-points are found to be un- 
stable and to coll apse into a Swe et-Parker current sheet in 
the MHD regime dBiskamrj, Il996t) . In a collisionless plasma 
X-point, stability can be maintained through coupling to a 



dispe rsive plasma mode dSturrockl 1966 



Shav and Drake 



1998). Recently it was discovered that X-points can be sta- 
bilized in the presence of MHD Hall effect so that the out- 
flow opejisjip OTi_Jarg^r_scale, thus making reconnection 



fast dShav e t al., 1998, 2004). On the other hand, most astro- 



physical plasmas are turbulent, and the heliosphere is most 
probably not an exception. This means that X-points can be 
created by turbulence and sustained at small scale by Hall- 
MHD effects, but only until turbulence itself collapses the 
X- points to form extended th i ck out flow regions (as observed 
by ICiaravella and Raymond! (12008b in multi-frequency ob- 
servations of solar flares). Even without turbulence it has 
been suggested that magnetic islands dynamically produced 
at X-points tend to be volume filli ng and to produce thick 
collisionless reconnection regions (lDrake et al.l |2006) with 
high reconnection speed. However, since turbulence is likely 
ubiquitous in astrophysics plasmas, we concentrate here on 
fast reconnection mechanisms in weakly stochastic plasmas. 

A model of fast magnetic reconnection that generalizes the 
Sweet-Parker scheme for the case of weakly stochastic mag- 
netic fields was proposed by ILazarian and Vishniac ( 1999 ) 
(henceforth LV99). Even though the notion of reconnection 
affected by turbulence is not new, in the LV99 model it is 
recognized that turbulence can decouple the width of plasma 



outflow region from the scale determined by Ohmic effects. 
In fact the outflow width is limited by the diffusion of mag- 
netic field lines, which depends on turbulence only (see cen- 
ter of Fig. [3), and can be much wider than the thickness of 
the individual current sheets (see bottom of Fig . [3j . Although 
reconnection events happen on small scales Am, where mag- 
netic field lines get into contact, a number of independent 
reconnection processes takes place simultaneously over ex- 
tended macroscopic current sheets L 3> An within a wide 
outflow region A ~ L. Therefore the effective reconnection 
rate is not limited by the speed of individual Sweet-Parker 
events on scale An (where plasma resistivity plays a domi- 
nant role), instead, it is enhanced by the large A, that depends 
on the magnetic field wandering. In such a situation it was 
found that reconnection speed is close to the turbulent veloc- 
ity in the plasma. In particular, assuming isotropically driven 
turbulence character ized by an injection scale I < L the re- 
conne ction speed is (ILazarian and Vishniaa 1 1999t ILazarian , 
2006) 



VrkVa 



1/2 



Yl 
v A 



(1) 



where Vi is the turbulent velocity at the largest scale and 
Va the Alfven velocity. Since turbulence in the heliotail is 
assumed to be weak, magnetic perturbations are compara- 
bly smaller with respect to the mean heliospheric magnetic 
field, therefore Vi < Va, i.e. turbulence in the heliotail is 
sub-Alfvenic. 

Numerical MHD calculations of wea kly stochastic mag- 
netic reconnection were performed by Kowal et aU ( 20091) 
and they proved that reconnection is fast and independent 
of Ohmic resistivity of the plasma, as shown in right panel 
of Fig |4] while it depends on the power of injected turbu- 
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lence (shown in left panel of Fig. @J. In these simulations, 
turbulence is preexisting and not related to reconnection pro- 
cesses. By varying the Ohmic and anomalous resistivity of 
the plasma, reconnection rate is not affected, confirming that 
in the presence of turbulence resistivity is not important. 

In the LV99 model, reconnection develops while the wide 
outflow region is filled with turbulent reconnected magnetic 
field lines movin g in opposite dire ctions. In fact, numeri- 
cal simulation by Kowal et al. ( 201 ll) show that the interface 
between the oppositely directly magnetic fields has a much 
more complex topology if compared to the Sweet-Parker re- 
connection mechanism, and also to the schematic representa- 
tion in Fig.[3](see Sec. |5]l. The outflow volume is filled with 
enhanced current density regions with loops of reconnected 
magnetic flux, where locally reconnection works faster since 
the current density reaches higher values. The magnetic 
loops shrink as a dynamical consequence of reconnection de- 
velopment, while multiple reconnection events happen at the 
same time due to the magnetic field stochasticity. 



5 Acceleration in reconnection regions 




v4 



Fig. 5. The simplest realization of magnetic reconnection at small 
scale and of acceleration as an energetic particle bounces back and 
forth between converging magnetic field lines. The converging ve- 
locity determines the reconnection speed Vr, while the advection 
of the accelerated particles entrained on the magnetic field lines, 
occurs at an outflow speed that in most cases is of order the Alfven 
velocity of the plasma Va- Particles bouncing at points A and B 
happens because either of streaming instability induced by ener- 
getic particles or magnetic turbulence in the reconnection region. 
In an actual turbulent plasma, the outflow region, at large scale, is 
filled with reconnecting loops and curren t sheets, each of which a 
possible acceleration site (see text). From|Lazarian (20051). 



Electric fields associated with reconnection events can ac- 
celerate energetic particles. For a particle of charge q, the 
typical energy gained in a reconnection process is of the or- 
der q(Vn/c)BAq, where X q is the coherence length of the 
particle within the reconnection layer. Efficient acceleration 
would require, therefore, both Vr, and X q to be large. How- 
ever, in general in any fast reconnection mechanism, the frac- 
tion of volume that is subject to resistive effects and reveals 
strong electric fields is small and most of the magnetic energy 
is converted into kinetic energy of the plasma instead. There- 
fore only a small fraction of the energy can be transferred 
through any fast reconnection process to energetic particles 
if direct electric field is involved. The observation of a large 
normal component of the electric field near an X reconnec- 
tion point in the Earth's magnetotail, was inte rpreted as Hall 



electri c field at the X-point current sheet by IWygant et al 
capable of accelerating ions to ~ 10's keV scale. 
In this paper we concentrate on the mechanism of energetic 
particle re-acceleration in weakly stochastic reconnection re- 
gions. 

In the LV99 mechanism, reconnection speed can approach 
Va, which can be appreciably large, therefore particles en- 
trained on reconnecting field lines bounce back and forth 
between the approaching magnetic walls while staying on 
the field lines that are contracting. This results in an in- 
crease of particle velocity with every bouncing, a s dis- 
cussed b y Ide Gouveia Dal Pino and Lazarian] ( 20031 12005 ) 
(see also iLazarianl ([2005)) where it was shown that recon- 
nection induces particle acceleration. The effect that indi- 
vidual magnetic loops shrinking in the reconnection region 
have on energetic particles, is equivalent to that of first order 
Fermi acceleration in magnetic mirrors. Fig. ^schematically 



represents the simplest realization of acceleration within the 
reconnection region expected within LV99 model. As en- 
ergetic particles bounce back and forth between converging 
magnetic fluxes, they gain energy. In the figure, the bouncing 
at points A and B is just an illustration of the process. In real- 
ity particles never pass by the same points in 3D, but they lo- 
cally stream along magnetic field lines and bounce back and 
forth through magnetic bottles that form in the reconnection 
region. 

The simple acceleration process represented in Fig. |5] 
can be easily quantified. An energetic particle with en- 
ergy E bouncing back and forth between a magnetic mir- 
ror will gain an energy AE ~ (Vr/c)E in every cycle. 
The process continues until particles gain enough high en- 
ergy to either diffuse perpendicularly out of the recon- 
nection region or get ejected by the outflow plasma at 
the Alfven velocity. This last possibility was c onsidered 
by Ide Gouveia Dal Pino and Lazarianl (120031 120051) . namely 
that particle diffusion perpendicular to the mean magnetic 
field is negligible. Perpendicular diffusion arises from mag- 
netic field wandering as particles scatter marginally perpen- 
dicular to the local magneti c field . Thi s effect was ac- 
counted by lYan and Lazarianl (2004, 2008) to describe cos- 
mic ray propagation, for instance. As mentioned, the proper- 
ties of turbulent plasma in the heliotail are not well know 
at this point, therefore it is difficult to quantify the dif- 
fusion regime, nevertheless, in general perpendicular scat- 
tering in sub-Alfvenic turbulence is fou nd to be subdom- 
inant with respect to p a rallel scattering dYan and Lazarian , 
1 2008 : Bere snvak et al.L 1201 II) (see also ILazarianl ([2006, 
2007)). In this case, the energy spectrum of accelerated 
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test particle is ( de Gouveia Dal Pino and Lazarian , 2003 



Lazarian and Opheri 2009) 

N{E)dE ~ E- 5/2 dE. 



(2) 



If perpendicular diffusion in the reconnection region is not 
negligible, particles bouncing between approaching recon- 
necting field lines of the magnetic mirror are not confined 
as within walls, but can keep bouncing while reconnection 
proceeds. In this situation particles may have cross field 
propagation but cannot escape from the large scale recon- 
nection region, producing a spectrum asymptotically reach- 
ing N(E)dE ~ E^dE. 

In case of re-acceleration of cosmic rays with a seed spec- 
trum E~ 2 - 7 , after acceleration it still becomes oc E~ 5 / 2 , i.e. 
harder than the initial spectrum. It is important to note that 
the expected energy spectrum Eq. [2] is an estimate based on 
a rather idealized situation. The derivation above considers 
only particles bouncing back and forth between the two re- 
connection layers. The actual picture of stochastic recon- 
nection in the LV99 model includes many simultaneous re- 
connection events happening at differe nt scales thr o ughou t 
the reconnection volume. Fig. l6l (from iKowal et al. ( 201 1 )) 
shows an evolved 2D snapshot of magnetic field configura- 
tion during reconnection from a nearly incompressible non- 
resistive MHD domain simulation without including kinetic 
effects, such as pressure anisotropy, the Hall term, or anoma- 
lous effects. The initial condition of the domain was set with 
eight Harris current sheets in a periodic box and a density 
profile corresponding to a uniform total pressure. A pertur- 
bation with random weak velocity fluctuations was used to 
enable spontaneous reconnection events. It is evident from 
the figure that the reconnection volume is filled with mag- 
netic loops (or islands) and that several reconnection events 
occur at the same times within the loops and along the current 
sheets between the loops. The simulation shows the exis- 
tence of merging loops with their resulting deformation and 
contraction, that provide appropriate conditions for particle 
acceler ation. This p i cture is very similar to the 2D simula- 
tion by iDrake etaD d2Q10h . where islands, or loops, are in 
fa ct only the 2D proj ections of 3D magnetic tubes, as shown 



m 



Kowalet alj (1201 lh 
The topological complexity of the reconnection region 
may have an influence on the actual spectral shape of the 
accelerated particles. Each local reconnection region (within 
a magnetic loop or a current sheet) provides the accelerated 
spectrum oc E~ b l 2 . But when energetic particles cross sev- 
eral reconnection regions they undergo further acceleration 
with a seed spectrum corresponding to that gained within 
the previously crossed reconnection region. It is well known 
from the theory of diffusive shock acceleration, that the sup- 
ply of a seed power law spectrum into an acceleration region 
leads to an amplification of the distribution without changing 



Fig. 6. Evolved 2D snapshot of magnetic field configuration where 
eight parallel Harris current sheets were perturbed in order to trig- 
ger plasma instabilities and reconnection events (represented in grey 
shades). 10,000 test particles, with initial thermal distribution with 
temperature corresponding to the sound speed of the MHD model, 
were injected in this plasma snapshot to study the acceleration 
mechanism induced by magnetic reconnection. The red and green 
colors correspond to regions where either parallel or perpendicu- 
lar acceleration occurs, respectively, while the yellow color shows 
locations where both types of acceleration occur. The parallel com- 
ponent increases in the contracting islands and in the current sheets 
as well, while the perpendicular component increases mostly in the 
regions between current sheets. The white boxes correspond to the 
sites where detailed d etermination of ac celeration properties were 
done (see text). From lKowal et alj J201 lh . 



the spectral index. However the increased number of high 
energy particles is accompanied by a decrease of the num- 
ber density at the low energy cutoff, leading t o a fl a ttenin g 
of the distribution at inte r mediate energy (see Bell (119781) : 
Melrose and Popd dl993l) : iGieseler and Jonesl (bood) ), and, 
therefore, to a harder energy spectrum. The exact slope of the 
final spectrum after bouncing within and escaping from sev- 
eral reconnection regions depends on the cutoff energy and 
on energy loss processes that particles undergo during ac- 
celeration and, most importantly, between acceleration pro- 
cesses. Scattering, for instance can degrade particle energy 
so that at each acceleration step a new low energy popula- 
tion is seeded into the process, leading to a softening of the 
spectrum. 

Non-linear effects from back-reaction of the acceler- 
ated particles might be important in reconnection pro- 
cesses. However so far the only approach to address back- 
reacti on of particles on re connecting plasma involved elec- 
trons (IDrake et all 120061) . In fact the evidence of back- 
reaction can be found in the si mulations of tes t particles 
propagating in the magnetotail dBirn et al. . |2004 and also 
in test particles studies in MHD models w ith magnetic is- 
lands dMatfhaeus et allll984tlKlierrIll994) . In lDrakeetaT 
(2006), where acceleration occurs in contracting loops 
formed in 2D collision-less reconnection, back-reaction is 
introduced by the term (1 — 87reji/i? 2 



where en is the ener- 



i.e. neglecti ng the back-reaction of accelerated particles, 
lLongain < ll992l) 



getic particle parallel energy averaged over the distribution of 
particle velocities. This would produce an accelerated spec- 
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Fig. 7. Top panel: kinetic energy evolution in time of 10,000 pro- 
tons in a 2D model of reconnection (i.e. with magnetic field on a 
plane). Bottom panel: the same distribution but in a fully 3D model 
of reconnection. In 2D the perpendicular component of particle ve- 
locities (in blue) becomes dominant over time with respect to the 
parallel component (in red). While in a 3D domain it is the parallel 
component to dominate. The energy is normalized to the proton rest 
mass. The background magnetized flow with multiple current sheet 
layers is at time 4.0 in Alfven time units in the model. Note that 
the transition from exponential energy grow to nearly linear occurs 
when the largest loop in the plas ma reaches t he siz e of a few tens 
the size of simulation box. From Kow al et"al] {201 ll) . 



trum ~ E 



3 / 2 , instead of the steeper 

Numerical simulations of test particles injected in the do- 
main represented in Fig. [6] with initial thermal distribution, 
show that particle velocity parallel (in red) and perpendicular 
(in green) to the mean magnetic field increases. Yellow color 
indicates the locations where acceleration increases both 
components without preference. But while in 2D domains 
perpendicular velocity mostly increases over longer integra- 
tion time, in 3D, where the loops develop in magnetic tubes, 
there is no such a limitati on making acceleration in the paral- 
lel dimension dominant dKowaletalll201lLl2012h as shown 
in Fig. [7] where velocities were sampled within the regions 
indicated with white boxes in Fig. [6] It is found that within 
contracting/deforming magnetic loops or current sheets par- 
ticles accelerate mainly through first order Fermi accelera- 
tion with pa rticles bouncing back and forth between converg 
ing mirrors (de Gouv eia Dal Pino and Lazarian I2003U2005 



Drake et all 1201(1 while outside these regions particles 
mostly undergo drift acceleration from magnetic field gra- 



dients. It is possible that turbulence far from loops, cur- 
rent sheets and diffusion regions, favors mostly second or- 
der Fermi acceleration mechanisms with particles being scat- 
te red by approaching and receding magnetic irregularities. 
In lKowal et alj ( 120121) it was argued that second order Fermi 
acceleration is the dominant process in purely turbulent plas- 
mas with no converging flow, although its rate is reduced. 
Moreover, reconnection layers in pure turbulence could be 
responsible for first-order Fermi acceleration of low energy 
particles. However, more studies are needed to fully un- 
derstand the interplay between different acceleration mecha- 
nisms in turbulent media. 

In shock acceleration, particles gain energy via plasma 
compression differences between the upstream and down- 
stream regions. Energy gain is described by Parker's trans- 
port equation jParkerl 1965 ), which was derived in the limit 
of strong scattering, and it is explicitly driven by plasma 
compression. On the other hand in magnetic reconnection 
it is possible to have acceleration even in an incompressible 
pl asma. In fact numerical MHD si mulations such as the one 
bv lKowal et al. (2009, 1201 ll 12012) were done in a nearly in- 
compressible regime, and as long as there is no strong scat- 
tering to maintain plasma isotropy the parallel energy gain 
dominates and, as a result, the particles entrained in the re- 
connecting magnetic flux gain net energy. 



5.1 Re-acceleration of Cosmic Rays 



The sectored heliospheric magnetic field regions in the helio- 
tail generated by the 1 1 year solar cycle is composed by 100- 
300 AU wide unipolar domains with turbulence scale likely 
of order 1 0- 1 00 AU, although the injection scale is not known 
precisely. The downstream solar wind motion in the heliotail 
induces converging flows in the turbulent plasma that ignite 
reconnection. As discussed in Sec. |4] turbulence creates the 
conditions for forming multiple simultaneously reconnecting 
magnetic fluxes (loops and current sheets) throughout large 
portions of the plasma. Acceleration takes place across the 
entire reconnection region and energetic particles are accel- 
erated through a sequence of multiple reconnection events 
(see Sec.|5j. Although second order Fermi acceleration from 
pure turbulence may occur as well, as long as reconnection 
is efficient first order Fermi acceleration is dominant. The 
overall process, therefore, takes place across regions that are 
comparable with the size of the unipolar magnetic domains, 
or even the size of the heliotail itself. The energy spectrum of 
those re-accelerated particles is, in the simplest case, ~ E 1 , 
with spectral index 7 = 5/2 or smaller as discussed in Sec. [5] 
which is flatter than the mean cosmic ray spectrum. Such ac- 
celeration can occur for as long as the cosmic ray particles 
are trapped within the reconnection regions. Using the sim- 
ple argument that the gyroradius should not be larger than 
the size of the magnetized region L zone , the maximum en- 



10 



Desiati & Lazarian: Cosmic Rays and Reconnection in the Heliotail 



ergy for a proton is (Longair 1992 1 



E,, 



0.5 



B \ { L z 



\ jiG i \ loo Ali) 



TeV. 



(3) 



The magnetic field strength in the heliotail is not known 
with pre cision, but we can assum e that it is of the order of 



1-4 fiG dPogorelov et al.L l2009al) . L zone is assumed to be 



within the range 100-500 AU (the higher bound being ap- 
proximately the heliotail thickness), therefore the maximum 
energy that cosmic rays can achieve is approximately be- 
tween 0.5 - 10 TeV. This means that the fractional excess 
region observed in the direction of the heliotail is likely ex- 
pected to have a harder spectrum than the average cosmic 
rays up to about 1-10 TeV. Above this energy the spectrum 
transitions back to the steeper ~ E~ 2J . Scattering processes 
within the heliotail can mitigate the acceleration effects and 
the related cosmic ray distribution at a given energy. Based 
on the global magnetic field structure in the heliotail, con- 
firmed in MHD simulations (see Sec. [3}, TeV cosmic parti- 
cles experience the lowest scattering along the line of sights 
parallel to the interstellar downstream flow, where magnetic 
field is strongly mixed at small scale. While away from this 
direction scattering in the unipolar magnetic domains scram- 
bles particles direction and effectively reduces the overall 
acceleration efficiency. At sub-GeV energies, the stronger 
scattering along the heliotail would degrade anisotropy an d 
spectral features. As stated in lLazarian and Desiatil (120101) . 
the properties of magnetic reconnections are still under ex- 
tensive study, and their level of complexity is being subject 
of debate. At the same time the very little explored tail region 
of the heliosphere, makes the problem under discussion here 
even more uncertain. However it is suggestive that the ob- 
servation of TeV cosmic ray arrival distribution and energy 
spectrum over small angular regions could be used to probe 
properties over the most remote regions of the heliosphere. 

Although there is no energy spectral determination in the 
sub-TeV energy range, the sig nificant hardenin g of the spec- 
trum obs erved by Milagro dAbdo et all 12008b and ARGO- 



YBJ dDi Sciascio and for the A RGO-YBJ Collaboration, 



20121) is indicative of a possible re-acceleration mechanism 
that involves a fraction of cosmic rays propagating from 
the direction of the heliotail. While waiting for other 
experimental results that can confirm a harder than average 
energy spectrum of cosmic rays within the localized excess 
region, the energy flux corresponding to the ^6 x 10~ 4 
fractional excess from ~ 10 GeV to a few tens of TeV can be 
estimated to be approximately between 10~ 9 and 10~ 8 erg 
cm~ 2 s , for 7 = 2.7 - 2.0, respectively. The corresponding 
average power dissipated in the re-acceleration of such 
energetic particles is approximately between 10 20 and 10 22 
erg s _1 . Even though a precise quantitative assessment of 
the power necessary to re-accelerate the fraction of energetic 
cosmic rays forming the Milagro localized excess region, is 
not possible at this point, this simple estimation shows that 



the fraction of heliospheric plasma power dissipated into 
cosmic ray kinetic energy is very small if compared to that 
transp orted by the solar wind (ss 10 27 erg s _1 , see iPark 
(fl%2k 

It is interesting to note that within the last few years exper- 
imental evidence that cosmic ray spectrum becomes harder 
at abou t 0.2-0.3 TeV/nuc leon has b een accumulat ed by 
ATIC-2 dWefel et alll2008h. CR EAM dAhnet all[20ich and 



PAMELA (Adr ianiet all 120111) . In particular the CREAM 



results seem to suggest that cosmic ray spectrum may be- 
come softer again at about 10 TeV/nucleon, although more 
observation is needed on this regard. The direct observation 
of a correlation between spectral features and arrival direc- 
tion would provide a breakthrough on the role of the heliotail 
in the TeV cosmic ray properties. 



6 Conclusions 

The observation that cosmic rays are anisotropic has gained 
special attention in the last decade, since it could provide in- 
formation about the galactic sources of the energetic particles 
and about the properties of the local interstellar medium and 
of the heliospheric magnetized plasma. Of particular interest 
is the evolution with energy of its angular structure, espe- 
cially of the tail-in anisotropy which appears as a broad ex- 
cess at sub-TeV energies from the direction of the heliotail, 
and seemingly degenerate into separate localized fractional 
excess regions above a few TeV. The directional coincidence 
of the tail-in excess at sub-TeV energies and of the most sig- 
nificant of the localized fractional excess regions at TeV en- 
ergies, with the heliotail provides a compelling connection to 
this little known extended portion of the heliosphere. 

Although we cannot exclude that other phenomena occur 
and might dominate the origin of the observation, such as the 
effect of energetic cosmic ray interaction with the turbulent 
ripples along the heliotail, in this paper another mechanism is 
discussed. Namely that a fraction of cosmic rays propagating 
through the heliotail are re-accelerated via first-order Fermi 
acceleration mechanism in weakly stochastic magnetic re- 
connection processes that originate in sectored magnetic field 
domains produced by the 1 1-year solar cycle. In general 3D 
numerical simulation show that such an acceleration mech- 
anism can be efficient up to a few TeV, where a flatter than 
average spectrum could arise, depending on the competing 
effects of multiple acceleration processes and escape or loss 
effects, and back-reaction. On the other hand the properties 
of magnetized plasma in the heliotail are not yet fully un- 
derstood, therefore details of cosmic ray propagation in this 
region are still uncertain. Sub-TeV cosmic rays may be ac- 
celerated over extended regions and may undergo more scat- 
tering, thus producing a broader arrival distribution. While 
multi-TeV cosmic rays undergo more efficient acceleration 
and their localized substructure in arrival direction are more 
related to the acceleration sites along the heliotail. Such ac- 
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celeration mechanism is intrinsically anisotropic and as long 
as scattering is sub-dominant it would generate a net energy 
gain that could explain the seemingly harder spectrum ob- 
served within the localized excess regions by Milagro. 

Acceleration processes in weakly stochas- 
tic magnetic reconnection regio ns as described 



by de Gouveia Dal Pino and Lazarian (2003) have been 
used in Lazarian and Ophei ( 20091) to explain the origin 
of the anomalous cosmic rays. The Voyager aircraft mea- 
surements show that the anomalous cosmic rays persist 
also downstream the termination shock, indicating that the 
site of their acceleration is within the heliosheath closer to 
the heliopause in the upstream interstellar flow direction. 
The sectored magnetic field arising from the 26 day solar 
rotation and originated by the difference between rotation 
and magnetic axes, are pushed away by the solar wind 
and compressed upstream toward the heliopause, causing 
magnetic reconnection and energetic particle acceleration. 
A similar model for the origin of a nomalous cosmic rays 
was proposed by iDrake et all ( 2010l) where the process of 
collisionless reconnection was discussed. In this paper we 
discussed a similar mechanism of the re-acceleration of en- 
ergetic cosmic ray particles, where the scale of the sectored 
magnetic field is significantly larger. The higher magnetic 
energy involved provides the possibility to accelerate higher 
energy particles in an observable manner in terms of a 
slightly anomalous energy spectrum and arrival distribution. 
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